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Abstract 

A four-node network consisting of a negative loop controlling a 
positive one is studied. It models some of the features of the p53 gene 
network. Using piecewise linear dynamics with thresholds, the allowed 
dynamical classes are fully characterized and coded. The biologically 
relevant situations are identified and conclusions drawn concerning the 
effectiveness of the p53 network as a tumour inhibitor mechanism. 



1 Introduction 

In the modeling of biological processes, dynamical networks play an impor- 
tant role. Metabolic processes of living beings are networks, the nodes being 
the substrates, which are linked together whenever they participate in the 
same biochemical reaction. Protein-protein as well as gene expression and 
regulation are also dynamical networks. 

A central role in the analysis of biological networks is played by the ori- 
ented circuits also called feedback loops [Snoussi & Thomas, 1993], [Thomas 
et al., 1995]. Feedback loops are positive or negative according to whether 
they have an even or odd number of negative interactions. The general un- 
derstanding is that positive loops generate multistability and negative loops 
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generate homeostasis, that is, the variables in the loop tend to middle range 
values, with or without oscillations [Snoussi, 1998], [Gouze, 1998]. 

An interesting situation, not much studied yet, arises when positive and 
negative loops interact and control each other. Here, one such network is 
studied. Besides its interest as a dynamical system, it also tries to capture 
some of the features of the p53 network. 

The p53 gene was one of the first tumour- suppressor genes to be identi- 
fied, its protein acting as an inhibitor of uncontrolled cell growth. The p53 
protein has been found not to be acting properly in most human cancers, 
due either to mutations in the gene or inactivation by viral proteins or in- 
hibiting interactions with other cell products (There are however a number 
of cases where apparently normal p53 does not to achieve tumour control). 
Apparently not required for normal growth and development, p53 is critical 
in the prevention of tumour development, contributing to DNA repair, in- 
hibiting angiogenesis and growth of abnormal or stressed cells [May & May, 
1999] [Vogelstein et al., 2000] [Woods & Vousden, 2001] [Taylor et al, 1999] 
[Vousden, 2000] . In addition to its beneficial anticancer activities it may also 
have some detrimental effects in human aging [Sharpless & DePinho, 2002] 

The p53 gene does not act by itself, but through a very complex network 
of interactions [Kohn, 1999]. The simplified network that is discussed in this 
paper, although not being accurate in biological detail, tends to capture the 
essential features of the p53 network as it is known today. In particular, 
several different products and biological mechanisms are lumped together 
into a single node when their function is identical. The network is depicted 
in Fig.l. The arrows and signs denote the excitatory or inhibitory action of 
each node on the others and the letters b,c,p,m,g denote their intensities 
(or concentrations). 

The p53 protein (p) is assumed to be produced at a fixed rate and to 
be degraded after ubiquitin labelling. The MDM2 protein being one of the 
main enzymes involved in ubiquitin labelling 1 , the inhibitory node is denoted 
(m). There is a positive feedback loop from p53 to MDM2, because the p53 
protein, binding to the regulatory region of the MDM2 gene, stimulates the 
transcription of this gene into mRNA. 

Under normal circumstances the network is "off" or operates at a low 
level. The main activation pathways are the detection of cell anomalies, 
like DNA damage, or abnormal growth signals, such as those resulting from 

1 MDM2 binding to p53 can also inhibit p53's transcription 
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Figure 1: A network with a controlled negative loop controlling a positive 
one: A simplified p53 network 

the expression of several oncogenes. They inhibit the degradation of the 
p53 protein, which may then reach a high level. There are several distinct 
activation pathways. For example, in some cases phosphorylation of the 
p53 protein blocks its interaction with MDM2 and in others it is a protein 
that binds to MDM2 and inhibits its action. However, the end result being a 
decrease in the MDM2 efficiency, they may both be described as an inhibitory 
input control (g) to the MDM2 node. 

The p53 protein controls cell growth and proliferation (c) , either by block- 
ing the cell division cycle, or activating apoptosis or inhibiting the blood- 
vessel formation (6) that is stimulated by several tumors and is essential for 
its growth. Therefore c and b form a positive loop. Expression of oncogenes 
or abnormalities is coded as an external control (g) that operates both on 
the p — m negative loop and the c — b positive loop. These two (controlled 
loops) will first be studied separately and later their interaction is taken into 
account. 

The dynamical evolution of the network is modeled by a piecewise-linear 
discrete-time system similar to the one used in [Volchenkov & Lima, 2003] 
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for other gene expression networks, namely 



p(t+l) = a p p(t) + W pm H(T m -m(t)) 

m(t + l) = a m m(t)+W mp H(p(t)-T p ) + W mg H(T g -g) 

c (t + 1) = a c c{t) + W cb H(b(t) - T b ) + W cp H(T p - p(t)) + W cg H(g - T g ) 

b(t + l) = a b b(t) + W bc H(c(t)-T c ) + W bp H(T p -p(t)) 

(1) 

H is the Heaviside function, the Tj's are the thresholds and in all cases 

a t + Y,W lk = l (2) 

k 

W ik > 0. This condition insures that all intensities remain in interval [0, 1]. 

A similar network was studied in [Vilela Mendes, 2003] modelled by a dif- 
ferential equation system. However, the piecewise-linear discrete-time model 
considered here allows for a more rigorous characterization of the dynamical 
features of the network. 

The study that we perform analyses the two basic network modules, 
namely the positive and negative loops. This type of modular analysis has 
been proposed before [Thieffry & Romero, 1999] as a tool to study biological 
networks. However this is an adequate methodology only if, in each case, the 
effect of the other modules is taken into account as an external control. This 
changes the usual conditions for multistationarity and stable periodicity. We 
will follow this extended modular approach, by considering the p — m loop 
as controlled by the g signal and the c — b loop as doubly controlled by g and 
p. 

2 The p-m system: a controlled negative loop 

p(t + l) = a pP (t) + W pm H(T m -m(t)) (3) 
m(t + l) = a m m(t) + W mp H(p(t)-T p )+W mg H(T g -g) 

The thresholds T p and T m divide the unit square [0, l] 2 into four distinct 
regions with different dynamical laws. This is illustrated in Fig. 3 (for the case 
g > T g ), where close to the corner of each region we have drawn the next-time 
map corresponding to the evolution of p and m. Because the Cli S clXG smaller 
than one, the dynamics tends to drive the system to a fixed point. However if 
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Figure 2: The p — m system. A controlled negative loop 



the fixed point of the dynamics is outside that region, the dynamical law will 
change whenever the system moves to another region. Therefore the system 
may never reach a fixed point and evolve forever. This is illustrated in the 
Fig. 4, where the regions and the corresponding fixed points are color coded 
for the two cases (g = or g — 1). Notice also that, for future reference, the 
four regions are coded as 0, 1, 2 and 3. 
Define the quantities 

fmp = Wmp (4) 

, W mg 
fmg = — 

Notice that from (j2J) it follows 

fmp + fmg = 1 (5) 

Examining Fig. 4 and considering all qualitatively different positions of the 
thresholds one obtains a complete classification of all possible dynamical 
behaviors of this system, which are listed in the following table: 
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Figure 3: The p — m system. The dynamical laws in each region for g > T g 
TABLE I. 
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Asympt. behavior 


Coding 


T m < fmg 
T m > fmp 


9<T g 


p — > 0, m — >■ fmg 


-»• 3 -> 2 

T 
1 


9>T g 


p — > 1 , m — > fmp 


2-^0^1 

T 
3 


T m < fmg 
T m < fmp 


9<T g 


p — >■ 0, m — > /rag 


-»• 3 -> 2 

T 
1 


9>T g 


oscillation 


0->l 

T 1 
2^3 


T m > fmg 
T m > fmp 


9<T g 


oscillation 


0->l 

T 1 
2^3 


9>T g 


p — > 1 , m — > /mp 


2 -»• -»• 1 

T 
3 


T m > fmg 
T m < fmp 


9<T g 


oscillation 


0->l 

T 1 
2^3 


9>T g 


oscillation 


0->l 

T 1 
2^3 



(6) 



Tp 




P 



Figure 4: Colour coded dynamical regions and the corresponding fixed points 
for the p — m system 

The period of the oscillations will depend on the specific values of the 
parameters, but the qualitative nature of the dynamics only depends on 
the relative positions of the thresholds and the quantities fmg and fmp. 
The coding represents the transitions between the regions defined by the 
thresholds. Depending on the numerical values of the parameters, the system 
may stay for a number of time steps in each region but then it is forced to 
follow the arrows in the coding scheme. 

Of the above four working regimens, the first two (T m < fmg) are the 
biologically relevant ones, because in the absence of oncogenes (g), normal 
p53 is expected to be at a low level. The other two regimens correspond to 
biologically detrimental action of this gene, provoking premature aging. 
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Figure 5: The c-b system: a doubly controlled positive loop 

3 The c-b system: a doubly controlled posi- 
tive loop 

b(t+l) = a b b(t) + W bc H(c(t)-T c ) + W bp H(T p -p(t)) (7) 
c(t + l) = a c c(t) + W cb H(b(t)-T b ) + W cp H(T p -p(t)) + W cg H(g-T g ) 

Here we consider this system as a positive loop, with dynamical variables 
b and c and two controls p and g. The relevant quantities are 

fcp = (8) 



c 



fob = 

fcg = 



L — a c 
L — a r . 



with 

fcp + fcb + fcg = l (9) 

and 

fbc = -p- (10) 
i — a b 

fbp = 

1 - a b 



with 

fbc + fbp=l (11) 

As before, the dynamics is characterized by analyzing the position of the 
fixed points for the 4 possible situations (p ^ T p , g ^ T g ). This is illustrated 
in Fig. 5, for particular values of the thresholds. The regions defined by the 
thresholds are, in this case, coded as a, b, c, d. We recall that the coding, 
listed in the tables, represents the required dynamical transitions between 
the regions. Depending on the numerical values of the parameters, the sys- 
tem may remain in each region a certain number of time steps before each 
transition. By adjusting the parameters we may therefore model different 
action delays in the network. The most important parameters controlling 
the delays are the dissipation parameters a«. 



Wcp Wcg 




Og<Tg P<Tp 
O g>Tg p<Tp 
g<Tg p>Tp 
<^ g>Tg p>Tp 



c 



Figure 6: Colour coded dynamical regions and the corresponding fixed points 
for the c — b system 
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TABLE 2 A. (Case g < T g and p < T p ) 



g <T g ,p<T p 


Asympt. behavior 


Coding 


T c < fcp 

J 6 < J°P 


6-> 1 

c > jcp -+- JCO 


a^d^b 

T 
c 


T c < fcp 
T b > fbp 


1 

c — > few + fcb 


a^b— >d 

T 


fcp + fcb >T C > fcp 
T b < fbp 


6-> 1 

c ^ fcp + fcb 


a— >c— >d 

T 
b 


fcp + fcb >T C > fcp 
T b > fbp 


b -> l,c -> fcp + fcb 
or 

6 -> /bp, c -> /cp 
or 

oscillation 


dO 
aO 
b^c 


T c > fcp + fcb 
T b < fbp 


6 -> /bp 

c^ fcp + fcb 


a^c^b 

T 
d 


T c > fcp + fcb 
T b > fbp 


b — > /bp 
c -> /cp 


b— >c— >a 

T 
d 



TABLE 2B. (Case g < T g and p > T p ) 



9 <T g ,p>T p 


Asympt. behavior 


Coding 




b -> 0,c -> 


aO 


T c < fcb 
T b < fbc 


or 

6 — > /be, c — > /cb 
or 


do 




oscillation 


b^c 


T c < fcb 


b^ 


c^b^a 

T 
d 


T b > fbc 




T c > fcb 


b^ 


b^c^a 

T 
d 


T b < fbc 




T c > fcb 


b^ 


c^a^b 


T b > fbc 




T 
d 
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TABLE 2C. (Case g > T g and p < T p ) 



g>T g ,p<T p 


Asympt. behavior 


Coding 


T c < fcg + fcp 
T h < fbp 


b^ 1 

c ->• 1 


c^d^b 

T 
a 


T c < fcg + fcp 
T b > fbp 


b^ 1 

1 


a^b— >d 

T 
c 


T c > fcg + fcp 
T b < fbp 


b^ 1 

1 


a^c^d 

T 
b 


T c > fcg + fcp 
T b > fbp 


&-> l,c-> 1 
or 

6 -> /ftp, c^ fcg + fcp 
or 

oscillation 


dO 
aO 
b^c 



TABLE 2D. (Case g > T g and p > T p ) 



g>T g ,p>T p 


Asympt. behavior 


Coding 


T c < fcg 
T b < fbc 


b^ fbc 

c-> fcg + fcb 


a^b^d 

T 
c 


T c < fcg 
T b > fbc 


b^fbc 
c -> fcg 


a^b^d 

T 
c 


fcg + fcb >T C > fcg 
T b < fbc 


b -> fbc, c ^ fcg + fcb 
or 

b -»• 0, c -> /cp 
or 

oscillation 


dO 
aO 
b^c 


fcg + fcb >T C > fcg 
T b > fbc 




c^b^a 

T 
d 


T c > fcg + fcb 
T b < fbc 




b^c^a 

T 
d 


T c > fcg + fcb 
T b > fbc 


6^0 


c^a^b 

T 
d 



Notice that the oscillations coded as b±=>c can only be realized if the dynamics 
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jumps directly between these two regions. Otherwise, if in the course of the 
dynamical evolution, the system crosses the regions a or d, it will be attracted 
to their fixed points. In this sense this oscillation code is not as robust as 
0^1 

the one ( t j. ) in the p — m system and might only be implemented for 
exceptional parameter values. 

4 Biological implications and the coupled sys- 
tem 

The analysis performed in Sections 2 and 3 is quite general and applies to 
any parameter region. However, not all parameter values will model biologi- 
cally normal conditions. For example, in the absence of cell abnormalities or 
oncogenes (g < T g ), normal p53 will be at a low level. Therefore, as follows 
from Table 1. 

T m < fmg (12) 

On the other hand, on the absence of oncogenes and low p53, that is, under 
non-pathological conditions, cell growth should not be explosive. Therefore 
this condition requires for the model 

T c >fcp + fcb 
T b > fbp 

as follows from Table 2A. 

With the above (normal condition) requirements we are now ready to 
analyze the behavior for the other regimens. First, it follows from (|13|) that 
also T c > fcb. Therefore one sees, from table 2B that if g < T g and p > T p , 
then b — > and c — > 0. That is, in the absence of oncogenes, expression of 
p53 completely blocks cell growth. 

On the other hand, one notices, from table 2C, that in all cases, when 
g > T g (presence of oncogenes) and p < T p , there is explosive cell growth 
(c — > 1). Only in the special case T c > fcg + fcp and T b > fbp there are 
other solutions, but even in this case c — > 1 is a possible solution. 

Finally, the most interesting situation is to analyze the effect of expressed 
p53 (p > T p ) in the presence of oncogenes (g > T g ) - Table 2D : 

- For T c > fcg + fob (and any T b ) p53 effectively controls cell growth, 
keeping it at a low level (c — > fcg). In this network model, a large value of T c 
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is related to the effectiveness of cell growth to stimulate the required increased 
level of blood supply. The larger T c the less effective is the stimulation signal. 

- For T c values smaller than fcg + fcb the inhibitory effect of p53 may not 
be so effective. For example, for fcg + fcb > T c > fcg and T b < fbc there 
are several solutions and the outcome will depend on the initial conditions. 
For T c < fcg and T b < fbc no control is possible. 

For parameter values of T m and T p for which the p — m system is driven 
towards a fixed point, the behavior of the coupled system may be read directly 
from the c — b tables of the previous section and the conclusions are as 
listed above. For oscillatory conditions on the p — m system, the behavior is 
more complex and several different regimens may be accessed. Here we will 
illustrate, by direct numerical simulation, the behavior of the coupled system, 
for regions of parameter values which, as discussed above, are biologically 
relevant (T m < frag, T c > fcp + fcb and T b > fbp). 




5 10 15 20 25 5 10 15 20 25 

t t 



Figure 7: Dynamical evolution of p, m, b, and c for fmg = 0.55, fbp = 
0.4, fcb = 0.35, fcp = 0.2, <n = 0.8, T p = 0.5, T rn = 0.5, T c = 0.6, T b = 
0.5, g > T g 
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Figure 8: Same as Fig. 7 with different initial conditions 

Figs. 7 and 8 show the dynamical evolution of the variables p, m, b, and 
c for g > T g and the parameter values 

fmg = 0.55, fbp = 0.4, fcb = 0.35, fcp = 0.2 

and 

T p = 0.5, T rn = 0.5, T c = 0.6, T b = 0.5 
This corresponds to the situation 

T m < fmg, T m > fmp 

in the p — m system and 

fcg + fcb >T C > fcg, T b < fbc 

in the c — b system. 

The parameter values are the same for both figures. One sees that de- 
pending on the initial conditions, in one case the p53 action keeps c at a low 
level and in the other it does not. 
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Figure 9: Same as Fig.7 with T m = 0.4 instead of T m = 0.5 



Figs. 9 and 10 use the same parameters for the c — b system, the only 
change being the value of T m that now is 0.4, that is, T m < fmg and T m < 
fmp. In this case the p — m system oscillates instead of converging to a fixed 
point, but the action on the c — b system is similar and it again depends on 
the initial conditions. 



5 Remarks and conclusions 

1 - A simplified four-node model displays biologically reasonable features 
which might, at least, serve as a toy model for the p53 network. The model al- 
lows for a complete mathematical characterization of its dynamical solutions 
and how they depend on the parameter values (thresholds and couplings). In 
biology, to change the threshold positions corresponds to stimulate or inhibit 
the mechanisms of gene expression. Therefore the dynamical identification 
of their role on the dynamics of the network might be a guide for therapeutic 
action. 
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Figure 10: Same as Fig.9 with different initial conditions 



2 - For very large values of the T c threshold, our model p53 (p) achieves 
control of cell growth in the presence of oncogenes. However, as illustrated 
before, there is a biologically reasonable region of parameters for which the 
successful action of p, as an inhibitor of cell growth, strongly depends on the 
initial conditions. In practice it means that it is only effective if the p level 
is already large enough before the positive feedback effect of c on b starts to 
play a role. 

It is known that in about half of human cancers p53 action is lost through 
mutation in the p53 gene. In many of the remaining tumors, the p53 gene 
is intact but it does not achieve the proper response [Woods & Vousden, 
2001]. Several mechanisms have been proposed to explain the loss of genet- 
ically normal p53 action. They include abnormal conformation , targeting 
for degradation by viral proteins, defective localization to the nucleus, am- 
plification of MDM2, loss of ability to inhibit MDM2 through proteins such 
as plA ARF or loss of kinases that phosphorylate MDM2 and p53. What our 
simple model suggests is that there is a range of biological parameters for 
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which success or failure of the p53 action is very much a question of chance, 
that is, it depends on the initial conditions of the process. This would be a 
purely dynamic explanation for the failure of normal p53. 

3 - The actual biological p53 network is highly complex and contains very 
many nodes and pathways [Kohn, 1999]. Our toy model has concentrated on 
extracting the abstract functional actions rather than reproducing the bio- 
logical detail. However, there are still many other dynamical mechanisms in 
the organic control of DNA abnormalities or pathological cell proliferation. 
A number of other genes either play a similar role or interact with p53. On 
the other hand some of the "defense" mechanisms of cancer cells are directly 
related to the operation of the basic network. For example, oxygen starving 
(hypoxia) of tumor cells by lack of blood supply (the b node intensity) acti- 
vates the CXCR4 gene and this activation causes the tumor cells to migrate 
to other organs (metastasis) [Staller et al., 2003]. Inclusion of this and other 
qualitatively different effects seems an interesting possibility. 

4 - Sections 2 and 3 contain a complete classification of the fixed versus 
periodic nature of the orbits. However, as mentioned before, the actual period 
of the oscillations depends strongly on the values of the parameters. In Fig. 11 
we illustrate this dependence by color coding the logarithm of the period of 
the orbits, in the negative loop, as a function of the parameters a m and a p . 
The parameters are chosen such that W pm = 1 — a m , W mp = 0.5(1 — a m ) and 
W mg = 0.5(1 — a m ). these conditions satisfy the oscillation condition and 
mantain normalization. 

To visualize the orbits corresponding to any parameter value a tool for 
Windows is avaiable for download at http : //www . ii . uam . es/ aguirre/p53_emu . zip . 

The behavior of the periods is rather discontinuous, switching 
from small to large values for small parameter changes. This occurs, 
inparticular , when some points of the orbit lie close to the thresholds. 

In Fig. 12 the rotation number i (the number of cycles divided 
by the period of the orbit, is depicted. This quantity displays 
a more continous behavior that the period of the orbits. Highest 
values of i correspond to high values of both a m and a p . However, 
there are still some discontinuities due mainly to orbits close to 
the intersection of the thresholds. 
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